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tions are presented and the validity of the First Law is investigated, and in these cases an expression 
for the energy is provided. 
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Qi I. INTRODUCTION 

U '. 

Recent observational data seem to imply an accelerating expansion of the visible universe leading to the 
so called Dark Energy issue. There are several possible descriptions of this acceleration. Among them, the 
. simplest one consists in introducing a small positive cosmological constant in the, up to now, unchallenged 

J> " framework of General Relativity (GR), the so called ACDM model. A generalization of this simple modifi- 

cation of GR consists in considering modified gravitational theories, where some combination of curvature 
invariants (the Riemann tensor, the Weyl tensor, the Ricci tensor and so on) replaces, or is added, to the 
classical Hilbert-Einstein action build on the Ricci scalar term R only. The simplest class of modified theories 
is i 7 '(i?)-modified gravity, in which the action is described by a function F(R) of the Ricci scalar R (for a 
00" review see Refs. \A 2]). Typically, these modified models admit de Sitter space as a solution, its stability 

having been investigated by several authors (see, e.g., [1-6]). Furthermore, viable F(R) models, namely 
those which are able to pass the local gravitational GR tests and also t o give rise both to inflation and to 
the dark energy epoch in a unified way, have been extensively discussed [7H10l|. Another interesting class of 
modified gravitational models, with roots in the principles themselves of gravitational theories, is given by 
Weyl's conformal gravity, where the quadratic Weyl scalar appears. We would like to stress that indeed the 
detailed study of Weyl gravity (as e.g. its role in Lovelock's theories [Til] , vacuum energy, dark matter, etc.) 
may definitely improve our understanding of these cosmological problems. 

One should note that the current acceleration of the universe is actually a low-curvature phenomenon, 
while it is usually recognized that black holes are associated with high-curvature effects. However, as stressed 
in Ref. [l2[ (see also the references therein) , the investigation of static spherically symmetric exact solutions 
in modified gravity models may have real interest also from the observational point of view with regard to 
strong gravitational lensing, where they could lead to some observable effects beyond those of GR, which 
could be used to prove or disprove these theories. 

In the present paper we focus our interest on the impact of these theories on the physics of the corresponding 
black holes (BH). Static, spherically symmetric (SSS) solutions describing static black holes in general 
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gravitational theories have been investigated in several papers. Typically, modified models admit de Sitter 
space as a solution, but the issue to find exact SSS metrics different from the Schwarzschild-dS one appears 
to be a formidable task since, already for very reasonable models, the ensuing equations of motion are much 
more complicated than the corresponding ones for GR. For that reason the number of exact non-trivial 
static black hole solutions known so far in modified theories of gravity is extremely small (see Refs. [l3l - 
UH). Furthermore, these black hole solutions are not expected to share the same laws of their Einsteinian 
counterparts. In parallel with ordinary BH solutions, among the physical quantities one should assign to 
these modified gravity black holes are a mass, horizon entropy, temperature, and so on. In particular, the 
issue of associating an energy (mass) to these black hole solutions is quite problematic (see Ref. [2pj : in M. 
Visser's terminology, most of the BH solutions in modified gravity theories are in fact "dirty black holes"), 
and the modern terminology usually refers to them as Lifshitz BH solutions. 

Several attempts addressed to find a satisfactory answer to the mass problem have appeared (see, for 
instance, Refs. [2l| - |23j . and references therein), but the definition of a sensible mass parameter is still a 



much debated question. In Ref. |24|, for a particular class of BH solutions, an identification for the BH 
energy has been proposed, as a quantity proportional to the unique constant of integration which appears 
in the explicit solutions of i 7 '(_R)-modified gravity models in vacuum. This identification was achieved by 
making use of a derivation of the First Law of black hole thermodynamics from the equations of motion 
of F(i?)-gravitv, by independently evaluating the entropy, via the Wald method (25[, and the Hawking 
temperature [26j, via quantum mechanical methods in curved space-times. Some explicit examples along 
these lines have already been discussed in Ref. [24| . 

Here we go further into this analysis by considering a particular class of (topological) SSS solutions in 



presence of Maxwell field [32j, and this inclusion is investigated in F(i?)-gravity and in Weyl conformal 



gravity models. With regard to this issue, we note that the Maxwell action and equations read, respectively, 

^ /•"'/•-., • V fc F**=0, (1) 



J-EM = 

where j,j,k, ... run from to 3, g is the determinant of the metric tensor gij, M. the space-time manifold, 
Vj is the covariant derivative operator associated with and the electromagnetic field strength. 

We will look for static, (pseudo)-spherically symmetric solutions (SSS) with various horizon topologies, 
thus the generic metric reads 

ds * = - e ^B(r)dt 2 + ^- +r 2 dal, dal = + P 2 d<t> 2 , (2) 

B(r) 1 — kp z 

where a(r) and B(r) are functions of r only. This form of the metric is chosen in order to deal with static 
spherically symmetric BH solutions. We recall that one has a BH solution if there exists a simple positive 
zero of B(r H ) = 0, with 2a ( r ») being finite and non vanishing. Dirty or Lifshitz BH solutions correspond 
to 2a( - rH) ^ 1. Where not necessary, the argument of all these functions will be dropped. Furthermore, one 
has a topological BH if the "horizon" manifold with metric da\ is a sphere S2, a torus T2, or a compact 
hyperbolic manifold Y2, according to whether k = 1,0,-1, respectively. 

Due to the spherical symmetry, it is easy to show that the non vanishing components of the electromagnetic 
field are 

P a i r )n P - a ( r )n 
F m =F tr = e —^-, F oi = _ c _Ji 7 (3 ) 

Q being the electric charge. As a result, within the SSS Ansatz, the Maxwell action simplifies to 

Iem = Yj dt J dr6 —^- ' ( 4 ) 

where Vu is the volume of the "horizon" manifold, namely Vi = Air (the sphere), Vq — |5t|, with r the 
Teichmuller parameter for the torus, and finally V-i — Ang, g > 2, for the compact hyperbolic manifold 
with genus g (see, for example [27j). 
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The paper is organized as follows. In Sect. II we obtain SSS solution in topological Maxwell-F(i?)-gravity. 
In Sect. Ill we extend to the electro- vacuum case the Weyl gravity topological SSS solutions of Refs. [28M30j | 
to observe that, also in this case, the solutions admit three free integration constants. In Sect. IV we discuss 
the existence conditions of BH event horizons. Sect. V is devoted to the evaluation of Wald's entropy in the 
well-know case of i 7 '(i?)-gravity and in the non trivial case of Weyl gravity. In this Section we also present 
a simple proof that the conformal invariance of Wald's entropy for conformally invariant gravity models. 
In Sect. VI, the topological version of the First Law of BH thermodynamics for F(.R)-gravity is presented, 
and a discussion regarding Weyl's conformal gravity is given. Conclusions are provided in Sect. VII. Finally, 
in the Appendix, a simple proof of the conformal invariance of the Killing static spherically symmetric 
BH temperature is presented. Throughout the paper we use units of k& = c = h = 1 and denote the 
gravitational constant Gm by n 2 = 8itGn, so that G^ 1 ^ 2 = Mpi, being Mpi — 1.2 x 10 19 GeV the Planck 
mass. Furthermore, we will also put 2k 2 = 1, for simplicity. 



II. TOPOLOGICAL SSS SOLUTIONS IN MAXWELL-F(R) GRAVITY 

The Maxwell-F(R)gravity model is described by the action I = Ir + I em, where the action for modified 
,F(i?)-theories reads 



I R=W 1 2 f d 4 xV=gF(R) =^ 1= f d A x^—g 

^ K J M J M 



F(R)-±F^F l3 



(5) 



(as advanced, we are useing units in which 2k 2 = 1). F(R) is an arbitrary smooth function of the Ricci 
scalar R which, using the metric in ([2]), assumes the form 

R' 9 R Oh 

R = —B" - — (4 + 3ra') - =y [(1 + ra' f + r 2 a"] + , (6) 

where here and in the following the prime ' represents derivation with respect to the r variable. 
The equations of motion for the action in ([5]) read 

( 9ij A - V,:V, + Ry) d ^±- - ~ 9ij F(R) = -l 9ij F rs F rs + ~ F ir F js g " . (7) 

Of course, only two of the latter field equations are independent since the metric depends on two arbitrary 
functions a(r) and B(r), only. One has 

F(R) = Rf-^L(k-B- rB') + f - (45 + rB') + 2Bf" + |J , (8) 

-(2/ + r/')-/" = 0, (9) 
r 

where, for the sake of simplicity, we have set 

/ - /(r) - dR ' f -dr~dr dR 2 ' and so on. 

In principle, given the function F(R), the above equations permit to derive both a(r) and B(r), that is, the 
explicit form of the metric. Alternatively, fixing the form of a(r) one derives B(r) and, as a consequence, 
the values of F and R as functions of r evaluated on the solution. From such expressions one can "try to 
reconstruct" the form of F as a function of R. In general, the Lagrangian F(R) one eventually finds is not 
unique, since one has to infer its form starting from the value it assumes on the solutions. For example, all 
Lagrangians of the form F(R) = Rf(R) with limfl^o Rf(R) = have the Schwarzschild solution. 

A comment is here in order. It has to be stressed that the general solutions of © and © are in fact 
physically acceptable only for some values of the parameters. In the following, even if not specified, it is 
always understood that all free parameters have been restricted to values which give rise to physical solutions. 
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Just to see how the reconstruction procedure works, we first take the derivative of equation © with respect 
to r, thus obtaining 

^ [rB'(4f - rf) + B(-4f + 2r 2 f")] - 2(a') 2 Bf + 

f f f" 20 2 

+ ^3 [4k -AB + 2B"r 2 } + ^ [2rB - 4B) + — [AB + rB'\ - = , (f 0) 

where, in order to simplify the result, we have also used Eqs. ([5]) and The usefulness of the latter 
equation is due to the fact that it only depends on / which, in principle, can be derived from (|9]), once a(r) 
is fixed. 

Here we shall consider two important cases, that is, a(r) =constant and a(r) ~ log(r 2 ) (of Lifshitz type, 
but in this case we will work with Q = 0). In the first case the metric reads 

dr 2 

ds 2 = ~B{r)di 2 + — + r 2 d ( r 2 1 (11) 

and equations (|6j)- (fT0|) notably simplify In particular, from © it follows that 

f"(r)=Q =*> f = ar + b, a,b = constant . (12) 

Putting this result into (fTU)) and integrating, one gets 

1) o = 0,6^0, B(r) = k+f+c 2 r 2 + £ 7 , 

2) a^0,b = 0, B (r) = ^ + ^ + c 2 r 2 + ^ (13) 

3) a, b ^ , B(r) = | + + c 2 r 2 + §^ + c ( £ - - - X- + r 2 log 

coi Ci, c 2 being arbitrary integration constants. The last expression is the more general form which the B(r) 
function can assume if the metric (|11[) is a solution of a generic F(R) theory with electromagnetic minimal 
coupling. We see that it depends, at most, on four independent parameters (a, b, cq, c 2 ). 
The corresponding scalar curvatures are given by 

1) a = 0,b ^0, R = -12c 2 , 

2) a^0,6 = 0, R=-12c 2 + ^-^, (14) 



3) a,6^0, i?=-12c 2 + 4+c -12 log 



ar+b j_ fc(12a J r- i +4ab 2 r+18a 2 br' ! ~f) J ) 



Now, in some cases we are able to "infer" the function F(R) which generates the solutions in (fl~3"|) . In the 
first case (a = 0, b ^ 0), using (fT3j) . (fT4|) and (jHJ one may recover GR 

a = 0,6^0, ^ R = 4A, (15) 

I F(i?) = (R - 2A) = , 

where we have set c 2 — ► —A/3 and 6 — > 1. This is the topological Reisnner- Nordstrom solution with non 
vanishing cosmological constant. We stress that the latter equation represents the value of F(R) evaluated 
on the solution. Of course, there are infinitely many functions F(R) which reduce (on shell) to F(R) = 2A 
for R = 4A and, in general, they do not give the correct solution. 
In the second case here considered, that is a 7^ 0, b — 0, we have 



a ^ 0,6 = 0. 



B(r) = l + ^+c 2 r 2 + £^ 

fl = _12 C2+ 4-f2;, (16) 

y F {R) = 2 f-§ pc , 

where a,ci,c 2 are arbitrary constants. In the special cases Q — or k = 0, namely in the absence on 
electromagnetism or in the fiat case, we are able to get r as a function of R. 
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If Q = 0, we get 



R + 12c 2 



F(R) = 2a y/k(R + 12c 2 ) • 



(17) 



The latter function satisfies the condition (TT21) and so we can try to use it in order to build up the action. 
One can directly verify, using ([7]), that the action 



-g d x 



M 



2a^/k(R+12c 2 ) 



has a vacuum solution of the form , with 



ci 



B{r) = 2 + ^ + C2r ■ 
In a similar way, if k = we have, with a < 0, 



-2Q 2 



5a(i?+ 12c 2 ) 



1/5 



Q 2 r-5 a (i?+12c 2 ) 

^ \R) = "rr 



2Q 2 



4/5 



(18) 



(19) 



(20) 



The latter function satisfies (|T2l and a direct computation using ([7]) shows that the corresponding action 



-g d x 



M 
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-5a(i?+ 12c 2 ) 
2Q 2 



4/5 



F %3 Fa. 



gives the required solution 



ci 



S(r) = -~ + c 2 r 



5ar 3 



(21) 



(22) 



In the third case here considered, that is a 7^ 0, 6 ^ 0, the method does not work owing to the presence of 
logarithmic terms which prevent the explicit reconstruction of the Lagrangian F(R), which formally reads 



F{R)= l - b R + ^ b+ t ar ^ 2 f 
y ' 2 2r 2 2a 2 r 2 (ar + b) 2 

However, if we choose cq = then the logarithm disappears from all equations and we get 



a^0,6^0, c = =• { R = -12C2 1 h 



fc(fe+4ar) 



(23) 



(24) 



F{R) = \bR^^ r 
These equations can be easily solved to get 

- * =^ F(R) = b(R + 6c 2 ) + 2a y/k{R + 12c 2 ) . 

R + 12c 2 

This is a topological charged generalization of the result found in Ref. [l!| . Also in this case one can directly 
verify that the field equations in (UJ are satisfied. 

Lifshitz like solutions 



Now we are going to consider the second case, that is a(r) ~ log(r z ), z being an (in principle) arbitrary 
parameter. With this special but important choice, Eq. ([9]) can be explicitly solved, to get 



}(z) = \Jz 2 + 10z + 1 , 



(25) 



G 



a,b being arbitrary integration constants. In the absence of an electric charge, that is Q = 0, and choosing 
6 = (or a = 0) also the Eq. (fT0|) can be explicitly solved, obtaining in this way a complicated expression 
for B(r), valid for any z. It depends on two arbitrary integration constants, say c± t 2, and for z ^ 1 it has 
the form 



70) 7± being given function of z, which read 

To 



a ^0,6 = 



a = 0,6 ^ 



7± = 

7o = 
7± = 



(z-1)Lp(2)-3(z + 1)1 
7»->M-i I 

4 ^ 

r 



p{z)-3{z + l) 



(z-i)1p(z)+3( z+i) ; 

7z + p(*)-l i 1 / z 3 -[pQ)+3] z 2 +[8p(z)+69] z-23p(z)+41 
4 1 !V p(z) + 3(z + l) 



(26) 



(27) 



(28) 



The corresponding metric can be considered as a generalization of the one of Clifton-Barrow (see below) . In 
the special case z = 1, a logarithmic term is also present. 

For example, for the following choices of the parameters, we get 



1) a ^0,6 = 0,z = l/2 




4k , _cj_ , 

7 T r 7/2 2 ' ; 

_9c 2 



F(i?) 



6fc 



8l£a 



(29) 



2) a^0,6 = 0,z = 2 




(30) 



3) a^0,6 = 0,z = l 



B(r) 
R 

F(R) 



Cir -(3+V3) _ 2fclog(r) 

3+v^ 3+V3 2 ' 

12fclogr+2(9+\/3)fc-6(3+73)c2 
(3+V3)r 2 

_ 4fe , (54+30V3)c 2 -12k log r 2+v/3 
T 3 " (S+v^r 3 



r 2+\/3 



(31) 



4) a = 0,6^0,z = l/2 



5) a = 0,6^0,z = 2 




rc 2 



if c 2 = , 
if k = , 



(32) 



(33) 



6) a = 0,6^0,z = l 



B(r) 
R 

F(R) 



-3+V3 



2fc log(r) 
-3+V3 



+ C 2 . 



-3+V3 

12fclogr+2(9-\/3)fc-6(3-v / 3)c 2 
(3-V3)r 2 
R _ 4k , (54-30v / 3)c 2 -12fclogr 
r r 3 ~T (3-\/3)r 3 



2- 



„2-%/3 



(34) 



The function F(R) for the two examples with z = 1 has not been explicitly written in terms of the Ricci 
scalar, because it is a quite complicated expression which we do not will use in what follows, while in the 
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example (4) we have also set C2 = (or k = 0), in order to satisfy Eq. ©. We observe that example (5) and 
the first part of example (4) are particular cases of the Clifton-Barrow metric, corresponding to 6 = 1/2 and 
6 = 1/4, respectively. 

In the previous section we have seen that a given metric can be the solution of different actions. In 
the same way, a given action can give rise to different SSS-metrics. For example, F(R) ~ R 3 ' 2 has the 
Schwarzschild solution, but also the solutions in the examples (4) and (5) above. 



Topological Clifton-Barrow solutions 



Here we present the topological generalization (without charge) of the Clifton-Barrow solution [13| . which 
is a SSS metric of the form considered in the previous section, with 



a(r) = lo 
B(r) 



5(1+2<S)/(1-<S) 



(1-Sf 



(I - 26 + 4S 2 ){1 - 25 - 28 2 ) 



1/2' 



{1 - 26 + 45 2 )(l - 26 - 25 2 ) 



C 



r (l-2<5+4<5 2 )/(l-<5) 



(35) 
(36) 



C being an integration constant and 6 ^ 1 an arbitrary parameter. 
The Ricci scalar reads 

_65(1 + £W1\ 

R - (26 2 + 26 -l) [y 2 ) ' (37) 



and, for 8 ^ —1, on the solution we obtain 

F{R) ~ R 1+s , (38) 



One can directly verify that (f36| arc solutions of the field equations corresponding to the Clifton-Barrow 
Lagrangian (|38|). The Clifton-Barrow metric is a particular case of the one considered in the previous section, 
Eq. (|2"6"1) . with one of the two constants (ci or c-i) dropped off and, in fact, it can be obtained using the 
method described above. 



III. CHARGED TOPOLOGICAL SSS SOLUTIONS IN WEYL CONFORMAL GRAVITY 

In this Section we generalize the higher gravity black hole solution of Riegert [28|, [2^] and its topological 
version [3(| in the Weyl model including the charge. We also present general solutions of pure Weyl gravity 
derived with the method of conformal transformations. 

To start, we write down the conformal invariant action of the model, in the form 

I = I W + Iem , (39) 
where the gravitational part is given by Weyl's conformal gravity 

I w =Zw I d 4 x y/^g C 2 , (40) 

Jm 

where w is a dimensionless parameter, usually restricted to the values w > 0, and C 2 — Cij rs C % i rs is the 
square of the Weyl tensor Cij rs . It is a scalar quantity which, in 4— dimensions, is given by 



C 2 = -R 2 - 2 H,, If + Ri jrs R ijr8 . (41) 
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The conformal gravity model has a very interesting future domain and its phenomenology has been investi- 
gated in Ref. [31|. 

Again, we look for SSS solutions of the form in ©. With this choice, one has C 2 = A(r) 2 /3r 4 , where 



A(r) = r 2 B" + (3m' - 2) rB' + 2 [rV + r 2 (a') 2 - rol + l] B -2k. 



The total action reads now 



I = V k J dt J dr e Q(r) 



wA(r) 2 
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2r 2 



(42) 



(43) 



We are dealing with a higher-order Lagrangian system, because the Lagrangian depends on the first and 
second derivatives of the unknown functions a(r) and B(r). The corresponding equations of motion (after 
simplification) read 



wA"B + wA' 



B B 1 



1 

4r 2 



w{A 2 + 4kA) 



0. 



(44) 



A" + A' 



0. 



(45) 



In order to find explicit solutions, we proceed as in the previous section, first considering the case 
a(r) ^constant. With this choice, one immediately gets 



A(r) = 6c + 



6&i 



and Eqs. (gH and (05) give the solution 
B(r) = k + 3c + — + ci r + c 2 r 2 , 



c (2fc + 3co) 



24w6i ' 



(46) 



CO) b±, C2 being arbitrary integration constants. The metric here obtained is the generalization to the charged 
topological case of the black hole Riegert solution [28|, H3| • 

The case b\ = has to be treated separately, because the limit in (|46j) is singular. In this case, we have 



B(r) 



Q 2 2 

h ci r + c 2 r . 

8w 



h = o, 



w 



±1. 



(47) 



ci, C2 being arbitrary integration constants. By conveniently choosing the integration constants of the Riegert 
solution one obtains the Schwarzschild, Schwarzschild-de Sitter(AdS) or the dc Sitter metric, correspondingly. 

In principle, one can try to find other solutions by solving the equations for a non constant a(r), but 
they are really complicated and so it is convenient to proceed in a different way, by using the fact that the 
Weyl-Maxwell action is conformally invariant. This implies that, if ds 2 is a solution of the field equations, 
then also ds 2 = il 2 (x k )ds 2 is a solution, ft being an arbitrary (smooth) function of the coordinates. In 
particular, starting from a SSS solution of the form @ and choosing Q — Q(r) > to depend on the radial 
coordinate only, one obtains again a static and spherically symmetric solution, which can be set in the form 
([2]), with the change of radial variable f = rO(r). In fact, 



ds 2 = ~Q 2 (r)e 2a ^B(r) dt 2 + + r 2 fi 2 (r) da 2 , 



and this assumes the form 



dS 2 



3 25(r) B{r) dt 2 



dr 2 

W) 



r 2 da 2 , 



(48) 



(49) 
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where, for simplicity, we have re-stated the original variable r by rfi(r) — >■ r, £l(r) —¥ r/H(r) and 



B(r) = B(E(r)) 



5(r) 
rS'(r) 



e 25 WB(r) = e 2Q ( H ( r »B(S(r)) 



1 2 



S(r) 



(50) 



In this way we have obtained a class of SSS solutions for the action (j40|) in the form ([2]), specified by the 
functions a(r), B(r). Any solution in this class is related to the original one by means of an arbitrary, positive 
(smooth) function H(r). As one can immediately see, S(r) — r corresponds to the unitary transformation 
f2(r) = 1< On the contrary, starting from Riegert's solution ([46)) and choosing an arbitrary S(r) ^ r one 
obtains a class of solutions with a(r) ^ given by means of the equation 



O 2o(r) _ 



S(r) 



[S'(r)]' 



dr |_S(r) 



(51) 



from which we see that, for fixed <S(r), one can in principle derive the function S(r) which provides the 
desired transformation. 



Riegert-Lifshitz solutions 

This is a class of solutions similar to the ones obtained for the F(R) case, but starting from that of Riegert, 
performing a conformal transformation as described in the previous section, and imposing a(r) to be of the 
form 



a(r) — log(7p) 



7 >0 : 



(52) 



where 7 is a dimensional constant and z an arbitrary parameter. 
From (|51[) , we immediately get 



S(r) : 
S(r) = - 



1 



7 log(gr) 



z = l, 



1 



q + jr z 1 j(z — 1) ' 



q being an integration constant. 

In the first case [z = 1) for B{r) we easily obtain 



B(r) 



Ci Ci 



log(qr) + (k + 3c ) log 2 (qr) - 761 log 3 (gr) , 



7 7 

while, in the second case (z ^ 1) and for simplicity choosing q = 0, for B(r) we get 



(53) 



B(r) = c + 61 h ci — + c 2 -5- , 



&1 - hPTjz ' 

- _ fc+3co 
C ° — (z-1) 2 ' 

ci- ^ 



(54) 



^ — C2 

c 2 — — j 



7(^-1) ' 



IV. BLACK HOLE SOLUTIONS AND EVENT HORIZONS 



The SSS metric may actually describe a black hole. We recall that an event horizon exists as soon as there 
exists a positive solution th of 



B(r H )=0, B'(r H )^0. 
The second condition describes a non extremal BH. 



(55) 
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For example, in the well known case of constant curvature, namely (fT5j) with Q = and k = or k = — 1, 
the conditions (1551) are satisfied if both c\ and the cosmological constant A are negative (AdS); on the 
contrary, if k = 1, one has a black hole solution whatever the sign of A if ci < again. A further restriction 
may be B'(th) > 0, namely the Killing surface gravity at the horizon to be positive. 

For the topological Clifton-Barrow solution one explicitly finds 

r H = { c/k)^ s y (1 - 2S+4S ^ , (56) 

which is positive if C / k > 0. In order to have a black hole, one has also to impose (1 — 26 + 4(5) /(l — S) > 0. 

Finally, in the case of the charged topological Riegert solution (|4"6")l , requiring positive surface gravity one 
has 

k + 3c + ci r H + c 2 r H + — = , c\ + 2c 2 th - > , r ff > . (57) 

If t>i < and C2 > 0, then it is easy to show that there always exists a positive root of B(rn) = 0, 
independently of the values of cq and c 2 ■ 

In the following we shall assume to deal with (topological) SSS black hole solutions satisfying conditions 
(|5"5f . with the additional requirement B'(rjf) > 0. 



V. THE WALD ENTROPY 



With regard to the evaluation of the entropy associated with a black holes solution, there exists a general 
method due to Wald [25[ , see also [13, H3| • The explicit calculation of the black hole entropy Sw is provided 
by the formula 



Sw = —2tt 



( SC 



r = r H 
t = const 



\5R 



i j rs 



\l2 dp d<p . 



(58) 



Here ^/h^dpdcj) — dVk represents the induced volume form on the bifurcate surface r = rjj ,t = const, while 
C = C(gij, Rijrs, ^kRijrs, ••■) is the Lagrangian density of an arbitrary theory of gravity. The antisymmetric 
variable ey is the binormal vector to the (bifurcate) horizon. It is normalized so that ey-e' 3 = —2 and, if the 
metric has the form in @, turns out to be 

e« = V-SdoflnW ~ $S}) = e tt « (<5°<5) - $5}) , 
Sj being the Kronecker symbol. As a consequence, one has 

SC 



(59) 



~>w 



2a{r H ) 



SR 



(60) 



Ah = VkT H being the area of the horizon. Note that the expression in is valid only because we are 
considering an SSS metric of the form ((2J. Furthermore, the variation of the Lagrangian density in (|58|) and 
(|60p has to be performed by considering Rij rs and gij as independent quantities. As an example, 



SR 



SR 



(g ir g JS 



9 zs 9 jr ) 



SR 



ijrs 



SR, 



0101 



,-2a(r) 



(61) 



For F(i?)-modified gravity the Lagrangian density has the form C — F(R) + Cem and, since the electro- 
magnetic Lagrangian does not depend on the Riemann tensor, one immediately obtains 



■>w 



A H f(r H ) 



f(r H ) 



dF(R) 



4 ' Jy ' 1 ' dR 

For pure Einstein gravity a — 0, f{rn) — 1 and one recovers the well known Area Law. 



(62) 
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A. Wald's entropy for conformally invariant gravity models 



In general the Wald entropy is not invariant with respect to conformal transformations, but when the 
action is conformally invariant, then also the entropy does not change. This can be easily seen by recalling 
that, under a conformal transformation as in (1481). one has 



n 2 



n 4 



(63) 



Rijrs — ^ 2 Rijrs + U (g , VO) . 



C = C(R 



tnnpq 3 



..) = n- 4 c, 



(64) 



where U is a function which does not depend on the Riemann tensor. If the action is invariant, both metrics 
gij and gij are black hole solutions and thus, for the metric in ([2]), we get 



Sw 



-2vr 



-2vr 



-2tt 



r = r H 
t = const 



r = r H 
t = const 



r = r H 
t = const 



V 9009119 



5C 



5R 



0101 



dp d<j> 



SC 



r=rn 

5C SRqioi 



^^0101 SRqioi 



dpd<j) 



5R> 



0101 



dp dcj) = Sw ■ 



(65) 



As it is well known, also the Killing surface gravity is a conformally invariant quantity and this means that, 
for example, all black holes described in Sect, (jllip have the same entropy and surface gravity. 

As an example, let us consider Weyl's conformal gravity. From (|40]) . one has C = 3wC 2 + Cem- Making 
use of dHJ and we get 



>V = -2AirwA H e 2a{rH) 
'I 



-24™ Ah 



2R awl - (g oa R n + g n R m ) + ± g^g^R 



, 2k 
3r / 3r 2 



(66) 



In the particular case of Riegert's solution (|46|) . the latter equation simplifies to 



Sw = — 487rwVfe ( — + c 



(67) 



VI. FIRST LAW FOR BLACK HOLE SOLUTIONS 

In Ref. an expression for the uncharged BH energy in a class of F{R) modified gravities has been 
proposed and identified with a quantity proportional to the unique constant of integration, which appears in 
the explicit solutions. The identification is achieved by making use of the derivation of the First Law of black 
hole thermodynamics from the equations of motion of i 7 '(i?)-gravity, evaluating independently the entropy 
via Wald's method and the Hawking temperature (2(| via quantum mechanical methods in curved space- 
times. In this section we will briefly revisit the case of F(i?)-gravity by discussing the results of Ref. [2J| to 
the topological Maxwell solutions. In the following, for the sake of simplicity we shall assume the charge Q 
as being a given quantity and not subject to thermodynamical variations. This, of course, is a restriction, 
implying that in the First Law the thermodynamical potential associated with the variation of the charge 
will be absent. 
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From the equation of motion ((5J) of F(i?)-gravity evaluated on the event horizon (namely B(th) = 0) and 
recalling the expression (|62p for the Wald entropy, we easily get 

B'ir H )^ = ^[(F-R f y H + 2kf] ■ Vk " 



It is well known [26| that all black holes have a characteristic temperature Tk related to the existence of 
an event horizon and, with the associated entropy, they have thermodynamical properties. For the static 
black hole described by the metric in ((2J, the Killing-Hawking temperature is proportional to the Killing 
surface gravity and has the explicit form 

e a ^B'(r H ) 

T K = ^ , ■ (69) 

This is a conformally invariant quantity (see for example [35L HE] • This make sense if the conformally related 
metric is a black hole solution with event horizon at r = rjj, and with regard to this issue a simple proof of 
conformal invariance of the BH temperature is presented in the Appendix. The above relation between surface 
gravity and BH temperature is a very robust and well known result (see for example [34j). The expression 
([55)) has been derived in several different ways, like the original Hawking derivation, by eliminating the 
conical singularity in the corresponding Euclidean metric or either by making use of tunneling methods 
recently introduced in Refs. [37], |38[ and discussed in details in several papers. 

As in general relativity, the temperature appears in a natural way in (|68p . Furthermore, if the entropy 
depends only on r# , as it happens in a large class of F(R) theories, Eq. (155)) can be used to derive the First 
Law of black holes thermodynamics and employed to define a specific BH energy [24[. Thus, we may write 



V k e a(r »*> 
16tt 



[(F - Rf)r 2 H + 2k f] dr H = T K dS w - pdV := dE K , (70) 



V=\v k r%, P-^^r, (*) 

where V is the volume of the black hole and p may be interpreted as the work term associated with the 
electromagnetic field. 

If the differential of the function in (|70|) is interpreted as variation due to an infinitesimal change of the 
size ru of the black hole, one may write 

Ek := ^ J dr H e^H) [ {F _ Rf y H + 2k f] |^ . (72) 

We recall that this is valid only if rn depends on a unique variable parameter, which will be proportional 
to the mass of the black hole. If th and, as a consequence, also the entropy depend on other variables, 
then other thermodynamics potentials will appear and the expression for the energy cannot be computed by 
means of the above expression. 

In what follows we shall apply these considerations to some exact solutions of F(i?)-gravity. We will 
sec that, in some interesting cases, Eq. (|72|) allows to identify the energy of the black hole with a suitable 
integration (mass) constant which appears in the explicit solution. 

As we have seen in (|13p . the function B(r) in general depends on three or four parameters, but in some 
cases some of them explicitly appear in the "reconstructed" Lagrangian and thus we can assume them to 
be a fixed constant of the model. Only the remaining parameters (the integrating constant of the solution), 
have to be varied in order to obtain thermodynamical equations. If only one free parameter survives, then 
the corresponding variation can be related to the energy via the first law of black holes, as explained above. 
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BH energy in F(R) gravity: Examples 



Here we explicitly compute the energy for some models we have derived in Sect|TIJ To start with , let us 
consider the first model in (|13l) . described by (|T5|) . Since F(R) = (R — 2 A) and Q is kept fixed, the only free 
parameter is c%, and from (IB"2"j) . (J5U), and ([72)1 . we get 

^=47' Tk = 4^ ' ^ = ^( fc -3 r ^ = 8^ Cl -32^' (?3) 

where we have used Bh = 0, namely 



*-£-r"« + 5r°- m 

Furthermore, for A = 0, one has fc = 1, and Ufe = 47r, £^ = Thus 
O 2 

ci=r H + ^-. (75) 
4r# 

This may be interpreted as the Christodoulou-Ruffini irreducible mass relation [39| , which allows to identify 
the BH mass M as M = 2c x 

As a second example let us consider the non trivial model in (|16p in the two special cases (Q = 0, fc 
and (Q ^ 0, k = 0). These are described by the actions in (|T8)) and in (|2ip . respectively, and the energy can 
be computed by the method described above. 

In the first case, entropy, temperature and energy, respectively, follows from (p2"j) . and ([721 . We 

obtain, 

4ff , m 1 / fc 4cA _ 3aVb 



S w = —akr H , T K = -—\ h — , E 1 ^ = — — ci . (76) 

47T 47T yrjj J 52-k 

We see that the entropy and energy are positive only if all parameters a, k, c\ have the same sign, but 
this means they must be all negative, in order to have a positive temperature, in particular k = — 1. On 
the contrary, the given constant ci has to be positive for having ru > 0, yielding an asymptotically AdS 
space-time. 

The second case in quite similar to the previous one. Now F(R) is given in (|20l) . and we get 
A H m 1 f-Q 2 4 Cl \ _ 3aV k V k Q 



■ clth: j Tk = — — a q- , Ek = c\ . (77) 

47r 47r \ ar H r 6 H J 16ir 32nrH 

Here we recall that a < 0. 

Asafurther example we consider the topological Clifton-Barrow solutions (|35 |) - (p6l) of model F(R) ~ i? (5+1 . 
One easily finds 

where we have introduced the dimensionless constant 

* 6 = / 2^-W(d-l)'(d + l)^ 1 \ 

Wl - 2J - 252 VI _ 26 + 4S 2 {25 2 + 2&-1)*)' 1 1 

In order to have the BH horizon and a positive Killing temperature, the range of the parameter S in the 
latter equation has to be restricted to the values already discussed. Some additional restrictions are also 
necessary to get a positive entropy. With such restrictions the energy in ((78)) is well defined and positive. 
As expected, in the limit <5 — >■ it reduces to the Misner-Sharp mass of general relativity. 
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Let us pause to summarize these noticeable results. In Eq. (|T0|) we have extended the validity of the 
Killing/Hawking temperature Tr-, which expression is derived in general relativity by using quantum me- 
chanics, to topological Maxwell- i 7 '(i?)-gravity. We stress that the Killing/Hawking temperature is well defined 
by the metric. In this way, we identify the black hole energy making use of the derivation of the first law of 
black hole thermodynamics from the equations of motion, by evaluating the related black hole entropy via 
Wald's method, and the Killing/Hawking temperature as the semiclassical result stemming from Hawking's 
radiation, as in the case of general relativity. 



A. Black hole Energy in Maxwell- Weyl conformal gravity 

We have seen in Sect. IIIII that the Lagrangian of Einstein- Weyl's gravity contains only a dimensionless 
parameter w, while the solution B{r) depends on three arbitrary integration constants cq, b\, c 2 - In this case 
the m coordinate of the black hole horizon will depend on several integration constants. In Ref. [24J, a 
particular case with Cq and c 2 held fixed has been considered. 

In the general case, when asymptotically the solution is anti de Sitter, the energy of the BH may be defined 
by means of the Euclidean action and the First Law of black holes thermodynamics holds with additional 
thermodynamical potentials (see Ref. for more details). With regard to this issue, let us imagine we 
are dealing with a BH solution, with th = Tu(c n ) being the larger positive solution of B{th) = 0, and c n 
constants of integration. Then, one has 

»=^=£(f^..+^»>§^..) ■* ~f|*.. m 

where the c n are all independent and B'{rn) is the partial derivative with respect to r evaluated on the 
horizon. Recall that Tk and Sw are computable, for example, from (|67p 

S w = -48irwV k (— + cq) . (81) 



Thus, we get 



T K dS w = 12T4 (ci dc - 4c 2 dh - 2b 1 dc 2 ) + dc , a = c °( 2fc + 3c °) . ( 82 ) 

Zwbi b\ 

In the particular case when Co and c 2 are held fixed, one recovers the result of reference [24| . namely 

E K = -48V k c 2 bx . (83) 

In the general case, all the constants of integration are true thermodynamic variables and, when it is possible, 
one has to make use of Euclidean methods in order to identify the energy. Then the First Law is shown to 
held true but with additional thermodynamic potentials, as explained in Ref. (40j . 



VII. CONCLUSIONS 



In this paper, we have considered several example of extended models of gravity coupled with the elec- 
tromagnetic field and we have studied charged topological static spherical black hole solutions, extending 
previous works. New black hole solutions have been presented. For modified F(R) gravity, we have inves- 
tigated the First Law of BH thermodynamics, and for a important class of topologically charged black hole 
solutions, and assuming the charge of the black hole as a given fixed quantity, we have obtained an F(R) 
BH energy formula. In the particular case of general relativity, our expression for the BH energy has been 
shown to be equivalent to a relation which defines the irreducible mass introduced by Christodoulou and 
Ruffini. 
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We have also considered the Maxell-Conformal gravity model. As we have shown in Sect. [V] the entropy 
is conformally invariant for Maxwell- Weyl's gravity and, since the temperature is conformally invariant too 
(see the Appendix), we have 

where the quantities on the right-hand side refer to the solution obtained by a conformal transformation 
of the kind considered in IIIII This means that, for example, all thermodynamical quantities related to the 
Riegert-Lifshitz solutions considered in Sect. IIIII can be computed directly by using the general formulas or, 
more easily, by replacing c„ with c„ coefficients by means of the relations in ([Mf . 

As a final remark, it is reasonable to check now in detail if the several examples of modified gravity models 
we have discussed in this paper do represent in fact feasible or viable models. For example, in the F(R) 
class, one should have [4l| 

dF(R) d 2 F(R) , x 

for R > Rq, Ro being the current Ricci scalar curvature. In the case of the Clifton solution, these inequalities 
are satisfied as soon as 6 is positive and small. The same check may be also performed for the other classes 
of F(R) models. 



Appendix 

In this Appendix we give a simple proof of the conformal invariance of the Killing-Hawking temperature 
for an SSS BH space-time of the form 

dg 2 = _ e 2«(r) 5(r)tft 2 + *i + r 2 ^2 (85) 

The event horizon is located at B(th) — 0, and the related Killing-Hawking temperature is Tk = e aH B' H . 
Let us perform a regular conformal transformation, with conformal factor Q (r), 

ds 2 = n 2 (r)ds 2 = ~e 2& ^B(r)dt 2 + ^— + r 2 da 2 , (86) 

B(r) 

where 

B(r) - ^rr , e 2 " (r) = ^e 2a ^ . (87) 

\V(r) 

In the new metric, since f = Q(r)r, the event horizon is still located at r = r#. As a consequence, the new 
Killing-Hawking temperature is Tk — e° H B' H . Now, 

B' 

B 'h = 752" ' ( 88 ) 

and 

e &H = n 2 H e aH , (89) 
and then it follows that 

f K = e aH B' H = T K . (90) 



Acknowledgments. LS has been supported in this work by a Short Visit Grant, to Barcelona, of the 
European Science Foundation on New Trends and Applications of the Casimir Effect. EE was supported in 



16 



part by MICINN (Spain), grant PR2011-0128 and project FIS2010-15640, by the CPAN Consolider Ingenio 
Project, and by AGAUR (Generalitat de Catalunya), contract 2009SGR-994, and his research was partly 
carried out while on leave at the Department of Physics and Astronomy, Dartmouth College, NH, USA. 



[1] S. Nojiri, S. D. Odintsov, Int. J. Geom. Me th. Mod. Phys., 4 , 115 (2007) arXiv: |hep-th/060l"213] S. Nojiri and 

S. D. Odintsov, Phys. Rept. 505, 59 (2011) [arXiv: 101 1.05441 [gr-qc]]. 
[2] T. P. Sotiriou, V. Faraoni, Rev. Mod. Phys., 82, 451 (2010). 

[3] G. Cognola, E. Elizalde, S. Nojiri, S. D. Odintsov and S Zerbini, JC AP 0502, 010 (2005). 

[4] G. Cognola and S. Zerbini, J. Phys. A 39, 6245 (2006 ) [arXiv:hep-th/0511233| . G. Cognola, L. Sebastiani and 

S. Zerbini, to appear in proceedings MG12. larXTv!T 006.1586 [ gr-qc] (2 009^ 
[5] V. Faraoni, Phys. Rev. D 74, 104017 (2006) [arXiv:astro-ph7 0610734 . 

[6] G. Cognola, M. Gastaldi and S. Zerbini, Int. J. Theor. Phys. 47, 898 (20 08) [arXiv:gr-qc/0701138 
[7] W. Hu and I. Sawicki, Phys. Rev. D 76, 064004 (2007) [arXiv:0705. 11581 [astro-ph]J. 
[8] E. Elizalde, S. D. Odintsov, L. Sebastiani and S. Zerbini, Eur. Phys. J. C 72 (2012) 1843 |arXiv:1108.6184l 
[gr-qc]]; E. Elizal de, S. Nojiri, S. D. Odintsov, L. Sebastiani and S. Zerbini, Phys. Rev. D 83 (2011) 086006 
|arXiv:1012 .2280 [hep-th]]; G. Cognola, E. Elizalde, S. Nojiri, S. D. Odintsov, L. Sebastiani and S. Zerbini, Phys. 
Rev. D 77, 046009 (2008) |arXiv:0712.4017l [hep-th]] ; G. Cognola, E. Elizalde, L. Sebastiani and S. Zerbini, Phys. 

Rev. D 83 (2011) 063003 [arXiv: 1007.4676 1 [hep-th]]. 

S. Nojiri and S. D. Odintsov, Phys. Lett. B 657, 238 (2007) [arXiv:0707.194T1 [h ep-th]]. 
S. Nojiri and S. D. Odintsov, Phys. Rev. D 77, 026007 (2008) [arXiv:0710.1738l [hep-th]]. 

D. Lovelock, J. Math. Phys. 12 , 498-501 (1971). 

S. E. Perez Bergliaffa and Y. E. C. de Oliveira Nunes, Phys. Rev. D 84, 084006 (2011) (arXiv: 1107.57271 [gr-qc]]. 
T. Clifton and J. D. Barrow, Phys. Rev. D 72, 103005 (2005) [gr-qc/0509059]. 
S. Deser, O. Sarioglu and B. Tekin, Gen. Rel. Grav. 40, 1 (2008) |arXiv:0705.166"9l [gr-qc]]. 
T. Multamaki and I. Vilja, Phys. Rev. D 74, 64022 (2006) [astro-ph /0606373 ; T. Multamaki and I. Vilja, Phys. 

Rev. D 76, 064021 (2007) [astro-ph/0612775| . 

S. Capozziello, A. Stabile and A. Troisi, Class. Quant. Grav. 25, 085004 (2008) arXiv:0709.0891 [gr-qc]]. 
T. R. P. Carames, E. R. Bezerra de Mello, Eur. Phys. J. C64, 113-121 (2009) arXiv:0901.0814 [gr-qc]]. 
R. Saffari and S. Rahvar, Phys. Rev. D 77, 104028 (2008) |arXiv:0708.T482l [astro-ph]]: Mod. Phys. Lett. A 24 

305 (2009) [arXiv:0710.5"635l [astro-ph]]. 

L. Sebastiani and S. Zerbini, Eur. Phys. J. C 71, 1591 (2011) arXiv:1012.5"230l [gr-qc]]. 
M. Visser, Phys. Rev. D 48, 5697 (1993). 

S. Deser and B. Tekin, Phys. Rev. D 67, 084009 (2003) hep-th/0212292 ; S. Deser and B. Tekin, Phys. Rev. D 



75, 084032 (2007) arXiv:gr-qc/0701140 



G. Abreu and M. Visser, Phys. Rev. Lett. 105, 041302 (2010). 

R. G. Cai, L. M. Cao, Y. P. Hu and N. Ohta, Phys. Rev. D 80, 104016 (2009) [arXiv:0910.2387l [hep-th]]. 

G. Cognola, O. Gorbunova, L. Sebastiani and S. Zerbini, Phys. Rev. D 84, 023515 (2011) |arXiv: 1104. 28141 

[gr-qc]]. 

R.M. Wald, Phys. Rev. D 48, 3427 (1993). 

S. W. Hawking, Nature 248 30 (1974); Commun. Math Phys. 43 199-220 (1975). 

L. Vanzo, Phys. Rev. D 56, 6475 (1997) [gr-qc/9705004] . 

R. J. Riegert, Phys. Rev. Lett. 53, 315-318 (1984). 

P. D. Mannheim and D. Kazanas, Astrophys. J. 342, 635-638 (1989). 

D. Klemm, Class. Quant. Grav. 15, 3195-3201 (1998) gr-qc/980805l]. 

P. D. Mannheim, Astrophys. J. 479, 659 (1997) astro-ph/9605085 ; Prog. Part. Nucl. Phys. 56, 340-445 (2006) 



astro-ph/0505266 . 

5. Habib Mazharimousavi, M. Halilsoy and T. Tahamtan, Eur. Phys. J. C 72, 1851 (2012) |arXiv:1110.5085l 

[gr-qc] ] • 

V. Faraoni, Entropy 12, 1246 (2010) arXiv:1005.2327 [gr-qc]]. 

M. Visser, Int. J. Mod. Phys. D 12, 649 (2003) [hep-th/0106lTl1 

T. Jacobson and G. Kang, Class. Quant. Grav. 10, L201 (1993) lgr-qc/9307002l. 

M. Nadalini, L. Vanzo and S. Zerbini, Phys. Rev. D 77, 024047 (2008) [arXiv:0710.2474l [hep-th]]. 
M. K. Parikh and F. Wilczek, Phys. Rev. Lett. 85, 5042 (2000). 

M. Angheben, M. Nadalini, L. Vanzo and S. Zerbini, JHEP 0505 , 014 (2005) [arXiv:hep-th/0503081j . 
M. Nadalini, L. Vanzo and S. Zerbini, J. Phys. A 39, 6601 (2006) [arXiv:hep-th/0511250| . ~ 



17 



[39] D. Christodoulou and R. Ruffini, Phys. Rev. D 4, 3552 (1971). 

[40] H. Lu, Y. Pang, C.N. Pope and J.F. Vazquez-Poritz, Phys. Re v. D 86, 044011 ( 2012) |arXiv: 1204. 10621 [hep-th] ] . 
[41] A. De Felice and S. Tsujikawa, Living Rev. Rel. 13, 3 (2010) [arXiv: 1002.49281 [gr-qc]]. 



